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Abstract. We explore the cosmological sensitivity to the amplitude of isocurvature 
fluctuations that would be caused by axions in the "anthropic window" where the 
axion decay constant /a ^ 10^^ GeV and the initial misalignment angle 0i <C 1. In 
a minimal ACDM cosmology extended with subdominant scale-invariant isocurvature 
fluctuations, existing data constrain the isocurvature fraction to a < 0.09 at 95% C.L. 
If no signal shows up, Planck can improve this constraint to 0.042 while an ultimate 
CMB probe limited only by cosmic variance in both temperature and iJ-polarisation 
can reach 0.017, about a factor of five better than the current limit. In the parameter 
space of /a and Hj (Hubble parameter during inflation) we identify a small region 
where axion detection remains within the reach of realistic cosmological probes. 
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1. Introduction 

The Peccei-Quinn (PQ) mechanism and concomitant axion is arguably the most 
plausible explanation for the smallness of the © parameter of QCD [1,2]. This case 
has become stronger in recent years because the experimental precision exploration of 
the quark mixing matrix reveals that the Kobayashi-Maskawa mechanism of explicit 
CP violation accounts perfectly well for all observations [3], so one can hardly argue 
that CP was a fundamental symmetry that is only spontaneously broken. Likewise, 
the possibility of a massless up-quark is now convincingly excluded [3], leaving few 
credible alternatives to the PQ mechanism. Perhaps the greatest advantage of the PQ 
mechanism is that it can be verified by the detection of axions. 

The CP- violating term in the QCD Lagrangian is of the form {as/STi)QG'^Ganyi 
with CKs the strong fine structure constant and G the colour field-strength tensor. With 
the identification © — > a//a the axion field a inherits its defining interaction with the 
axial QCD anomaly. The axion decay constant /a is the primary parameter for axion 
physics. In particular, all couplings are inversely proportional to /a and the axion mass 
is determined by 

mj^ = ^"^""^'^ m^f^ = Aa. (1.1) 
mu + md 

Here = 135.5 McV is the pion mass, /^r = 93 MeV its decay constant, and the 
up- and down-quark masses with Z = m^/m^ = 0.35-0.6 [3]. A series of overlapping 
experimental and astrophysical constraints suggests /a > 10^ GeV [3,4], so axions are 
very light and "invisible," in any case far more elusive than neutrinos. 

However, during the QCD epoch of the early universe, a non-thermal mechanism 
produces axions as nonrelativistic coherent field oscillations that can play the role of 
cold dark matter [5]. In terms of the initial "misalignment angle" ©i = ai//a relative to 
the CP-conserving minimum of the axion potential, the cosmic axion density is [6] 

^a/^^- 0-195 ©?(^^j . (1.2) 

If ©f is of order unity, axions provide the dark matter of the universe if /a ~ 10^^ GeV 
(ma ~ 10 /leV). As such they can be detected by experiments that exploit the generic 
axion-two photon interaction, allowing for the coherent conversion into excitations of 
a high-quality microwave cavity placed in a strong S-field [2,7]. If axions are in the 
mass range 1-100 /leV and form the dark matter of our galaxy, the ADMX and New 
CARRACK experiments are poised to find them within a decade. Such a discovery 
would both reveal the nature of dark matter and verify the PQ mechanism for CP 
symmetry restoration. J 

Meanwhile one may study other "windows of opportunity" where axions could leave 
a detectable trace. Ongoing efforts include hehoscope searches that are sensitive to solar 

I In view of the fundamental importance of such a measurement it is perplexing how little worldwide 
effort is directed towards this goal compared with the vast range of WIMP dark matter searches. 
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axions if /a ~ 10^ GeV, corresponding to sub-eV masses [8, 9]. Of course, if solar axions 
were found something would have to be wrong with the supernova 1987A energy-loss 
limit that requires /a > 10^ GeV. For /a < 10® GeV axions thermalise in the early 
universe after the QCD epoch, so sub-eV mass axions would provide a hot dark matter 
component [10-12]. If future surveys were to reveal a hot dark matter fraction beyond 
the minimal neutrino component, an axion interpretation would be conceivable. 

We consider here axions in another range beyond the classical cosmological window. 
In a scenario where the PQ symmetry is not restored during or after inflation, a single 
value — TT < ©i < -l-TT determines the axion density. It is possible that ©i <^ 1, allowing 
for /a ^ 10^^ GeV. This case is motivated because the PQ mechanism presumably is 
embedded in a greater framework. In particular, the PQ symmetry emerges naturally in 
many string scenarios, but then it is hard to obtain /a much below 10^^ GeV [13]. The 
axion is the Nambu-Goldstone boson of the global f/(l)pQ symmetry that is broken at 
a scale wpg = N fg_ where N is an integer. (Both a/upQ and a//a must be 27r-periodic 
because a/vpq is the phase of a new Higgs field responsible for spontaneous symmetry 
breaking and the QCD vacuum is periodic in © = a/ fa-) Since — tt < a/vpq < -|-7r 
encounters N minima of the axion potential, there is a formidable cosmological domain- 
wall problem unless = 1 or inflation occurs after PQ symmetry breaking, further 
motivating a "late-inflation scenario" [5]. 

Postulating ©i ^ 1 may look like shuffling the QCD problem of a small © parameter 
into an unnatural initial condition for our universe. However, a GUT-scale value of 
/a ~ lO^'' GcV would only require a modest fine-tuning of ©i < 3 x 10^^ [14]. More 
importantly, anthropic selection suggests that the cosmic baryon/dark matter ratio can 
not be too small, so for any /a the most probable value for ©i is near the one required 
for the observed dark- matter density [15, 16]. 

Following these papers we emphasise that in our late- inflation scenario ©i, and 
therefore the axion dark matter density, is an environmental parameter of our universe, 
not a fundamental one, and that its prior probability distribution is known to be flat on 
the interval — tt < ©i < +tt. This contrasts with typical other cases of anthropic 
reasoning where usually it is not known if a given cosmological or particle-physics 
parameter is fundamental or environmental. In the axion case anthropic reasoning 
is unavoidable to quantify if a given value of ©i is natural or not. 

Therefore, the "anthropic axion window" with ^ 10^^ GeV is well motivated. 
In this case ma <^ 10 /leV and axion dark matter is difficult to detect with the cavity 
technique because of the large required cavity size and magnetic field region. § An 
alternative signature is provided by primordial isocurvature fiuctuations that can show 
up in future data. In our late-inflation scenario, the massless axion held is present 
during inflation and thus acquires the usual amplitude fluctuations imprinted by the 
de Sitter expansion. When axions acquire a mass during the QCD epoch, these 
fluctuations become dynamically relevant in the form of isocurvature fluctuations that 

§ Ideas for overcoming this limitation were discussed by S. Thomas in a talk given at the workshop 
"Axions at the Institute for Advanced Study" (20-22 Oct 2006, IAS, Princeton, New Jersey). 
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are uncorrelated with the adiabatic fluctuations inherited by all other matter and 
radiation from the inflaton flcld. The isocurvature amplitude depends on both /a 
and Hi^ the Hubble parameter during inflation, so observational limits on isocurvature 
fluctuations exclude certain regions in this parameter space [17-22]. 

Since there is no trace of isocurvature fluctuations in existing data, perhaps a more 
interesting question is the remaining window for axions to show up in future. One 
forecast was recently provided by the CMBPol Study Team Collaboration [23]. The 
purpose of our paper is to provide a more detailed sensitivity forecast on the amplitude 
of isocurvature fluctuations. We explore possible degeneracies that may exist between 
the isocurvature fluctuations and other cosmological parameters in both present and 
future data, and determine the region in the underlying axion parameter space that is 
realistically within the reach of future cosmological probes. 

To this end we review in Sec. 2 the axion dark matter and isocurvature contributions 
in terms of underlying physical parameters. In Sec. 3 we derive our constraints and 
sensitivity forecasts on the axion-type isocurvature fraction. We interpret these results 
in terms of underlying physical parameters in Sec. 4 and conclude in Sec. 5. 



2. Cosmological imprint of axions 

2.1. Energy density 

To specify the cosmological imprint of axions we begin with their contribution to the 
matter density of the present-day universe. When the PQ symmetry breaks at some large 
temperature T ~ fpq, the relevant Higgs fleld will settle in a minimum corresponding to 
©i = c^iZ/ai where — tt < ©i < +7r. We assume that this happens before cosmic inflation, 
so throughout our observable universe we have the same initial condition except for 
fluctuations imprinted by inflation itself. When the universe cools and evolves towards 
the QCD epoch, instanton effects become important and produce a potential for the 
axion field. The corresponding temperature-dependent axion mass is [6] 

m,{T) ~ I ^i-t (g^) ' 1^- ^otT> Aqcd, (2.I) 
I Aa//a for T < Aqcd, 

where Anst - 3.964 MeV^ n ~ 6.878, Aqcd ^ 380 MeV, and A^ ~ 77 MeV, assuming 
Z = 0.5. 

The axion field begins oscillating at 3if(Tosc) — ^a(^osc)- The Hubble parameter 
during radiation domination is 

H''(T)^^^—gJT)T\ (2.2) 

where Mpi = 1.221 x 10^^ GeV is the Planck mass. For T > Aqcd the effective number 
of thermally excited degrees of freedom is 5'*(Tosc) — 61.75, whereas g*{Tosc) — 10.75 for 
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T <S Aqcd- Together equations (2.1) and (2.2) give 

f 9.16 X 10^ MeV (tot^^^)""'''' for T^^^ > Aqcd, 

-tosc — S , / # \ -0.5 l^-<jj 

[ 6.65 X 10^ MeV [joJ^) for Tosc < Aqcd- 

We see that the hmit T^sc ^ Aqcd is only relevant for /a ^ 3 x 10^^ GcV. We will 
usually focus on /a values around or below 10^^ GeV and therefore concentrate on the 
case where the axion field begins oscillating during the QCD epoch. 

The traditional calculation of the present-day axion density [5] was recently 
revisited in detail, taking into account modern values for the relevant QCD parameters 
and anharmonic corrections for the axion potential [6]. However, since we are interested 
in the anthropic axion window where the initial misalignment angle is small, anharmonic 
effects are not important. We also ignore the possibility of axion dilution by late entropy 
production, so 

a;a ^l^a/^^^ 0.195 (^^^-^j (6,^ + 4). (2.4) 

Here, the initial misalignment angle ©i is to be interpreted as the average in our Hubble 
volume, ©i = (©)i, while 

4 = 4j 

is the inflation-induced variance, with Hi the Hubble parameter during inflation. In 
other words, the relevant initial parameter is (0^)i = ©f + o'q. 

We will usually assume that all of the cold dark matter consists of axions, so 
according to current cosmological data cUa = — 0.11 [20]. Assuming as is small, we 
then find 

©i ~ 0.75 ( J j (2.6) 



as a unique relationship between the initial misalignment angle and the axion decay 
constant. 



2.2. Isocurvature fraction 

Since the axion is essentially massless during infiation, quantum fluctuations are 
imprinted on the axion field in the same manner that the inflaton fleld acquires its 
perturbations. Because of the axion's negligible contribution to the total energy budget 
of the universe during inflation, these fluctuations do not perturb the total energy density 
and manifest themselves only as perturbations in the ratio of axion number density to 
entropy density, i.e., 5p — 0, and 5{ng_/ s) ^ 0. For this reason they are known as 
entropy or isocurvature perturbations. As the axion field provides some or all of the 
cold dark matter of the universe, these perturbations contribute to the temperature and 
polarisation fiuctuations in the CMB. 
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Axion- induced isocurvature fluctuations are uncorrelated with the adiabatic 
fluctuations inherited by other matter and radiation components from the inflaton. 
The spectrum of entropy perturbations is 

m - ^^j^- (2-7) 

Assuming 60 = Q — (0) has a Gaussian distribution, the perturbation power spectrum 
can be written as 



2al (20f 



2 • 



where 6i and ctq are the mean and the variance of the initial O distribution. As we 
shall see later, for all of the viable parameter space, 0f ^ (7@, thus 

In the second equality we have approximated the power spectrum by a power law with 
spectral index Ui^o = 1 — 2e ~ 1, where e is the first infiationary slow- row roll parameter 
and fco = 0.002 Mpc^^ the pivot scale. Likewise, adiabatic perturbations sourced by the 
inflaton fleld are encoded in the curvature power spectrum, 

77-2 f U\ "ad-l 

where is the adiabatic spectral index. 

Because fluctuations of the axion fleld are uncorrelated with those of the inflaton 
fleld, the primordial scalar power spectrum is simply the incoherent sum of the adiabatic 
and isocurvature perturbation spectra, 

v{k) ^ {\n{kf\) + {\s{kf\) . (2.11) 

The resulting CMB temperature and polarisation anisotropy spectra are given by a 
similar sum, but with the appropriate transfer functions Tx{k) folded in [X — T,E), 

dk ^„ J , , , ^„ J / ,n,i„/,noh f dk 



oc 



/ y T^^k)Tf{k){\n{kr\) +J- Tr{k)T-^{k){\S{kn. 

(2.12) 



for 02 > a2 



Here it is useful to define the isocurvature fraction a 

- {m-n + {\sikn -\^^ for 0? « ai, ^ ■ ^ 

where = V{k = ko) is the amplitude of the total primordial scalar power spectrum 
at the pivot scale k^, and we have assumed a <^ 1. Given 10~^, the solution 

0f -C ctq cannot be realised. In the opposite limit Qf ^ cTq we find 

(2.14) 

where we have used equation (2.4). 
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2.3. Gravity waves 

The misalignment mechanism per se does not produce gravity waves. However, all 
inflationary scenarios produce some amount of tensor perturbations, which is in principle 
manifest in the CMB temperature and polarisation anisotropics. The quantity of interest 
is the tensor-to-scalar ratio r = At/^Sj defined as the ratio of the primordial tensor 
to the scalar perturbation spectrum sX k = k^. The amplitude of tensor perturbations 
from inflation is A^: = 16i/|/(7rM|j). Thus, the Hubble expansion rate during inflation 

Hi^ ^^TrAsr Mpi 

~ 1.33 X 10^^ GeV ( , , ^ , J (-— ) (2.15) 
is directly probed by a measurement of r. 

3. Present constraints and future sensitivities 

3.1. Observational signatures of axion-type isocurvature perturbations 

To discuss the observational imprint in the CMB we recall that the axion field causes 
scale-invariant CDM isocurvature fiuctuations that are uncorrelated with the dominant 
infiaton-induced adiabatic fiuctuations. Compared to adiabatic initial conditions, the 
growth of isocurvature CDM density perturbations with wavenumber k is suppressed 
by a factor of k"^ on sub-horizon scales [24]. Hence, if we are looking for traces of a 
subdominant isocurvature perturbation with a similar primordial scale dependence as 
the dominant adiabatic mode, the signal will be strongest at large scales. We illustrate 
this point in Fig. 1, where we plot the CMB angular power spectra for a model with 
purely adiabatic perturbations {a — 0) and riad = 0.96, and for a corresponding model 
with purely CDM isocurvature perturbations (a = 1) with riiso — 1- For uncorrelated 
isocurvature perturbations, cases with < a < 1 simply correspond to weighted sums 
of these two spectra, so even for largish values of a, the spectrum at multipoles i > 200 
will be completely dominated by the adiabatic contribution. 

It is also useful to contrast the isocurvature signal with the cosmic variance, which 
poses a fundamental limit to the amount of information to be extracted from CMB 
anisotropics. In Fig. 2 we compare the contribution AC^ of a subdominant isocurvature 
part to the total Ci with the intrinsic uncertainty due to cosmic variance, Agy = 
Cf-^Cj^ -sj , with X,Y E [T, E] . This shows that data of small-scale perturbations 
by themselves carry hardly any information about axion-type isocurvature. 

Based on existing measurements of the CMB anisotropics, the main signature of 
scale-independent uncorrelated CDM isocurvature, i.e., an enhancement of the Sachs- 
Wolfe plateau, can be mimicked to a certain extent by changing the values of other 
cosmological parameters. In particular, one can expect a to be degenerate with the 
total matter density fl^ and the adiabatic spectral index nad- Adding data of the large- 
scale matter power spectrum and of distance probes such as supernovae will however 
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i 

Figure 1. CMB angular power spectra for purely adiabatic initial conditions (a = 0, 
dark red lines) with riad = 0.96 and for purely CDM isocurvature initial conditions 
(a = 1, light red lines) with a scale invariant spectrum riiso — 1- All other cosmological 
parameters are kept fixed. 




10 100 1000 2000 

•C- 



Figure 2. Scale dependence of the "signal to noise" for isocurvature fractions of 
ct = 0.1 (blue lines), and a = 0.001 (pale blue lines). Plotted is the distortion ACe 
of the CMB angular power spectra caused by adding a subdominant scale-invariant 
uncorrelated CDM isocurvature component, given in units of the cosmic variance Acv- 
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allow us to break these degeneracies and thus improve constraints on a indirectly. 



3.2. Constraints and forecast 

For inferring constraints from present data and predicting the discovery potential of 
future CMB experiments, we consider a cosmological model with seven free parameters: 
the isocurvature fraction a plus the usual six parameters of the "vanilla" ACDM model. 
These parameters and their benchmark values to generate mock data for our forecasts 
are the baryon density c^b = 0.0226, the dark matter density Uc = 0.109, the Hubble 
parameter h = 0.715, the redshift to reionisation z^e = 11.36, the scalar spectral index 
Had — 0.962, the normalisation ln(10^°74s) = 3.195, and the isocurvature fraction a — 0. 
The spectral index of the primordial isocurvature power spectrum is fixed to riiso = 1- 
For the other parameters we adopt flat priors in the data analysis. 

We do not vary the tensor-to-scalar ratio r, since in the anthropic parameter region 
any possibly observable amount of tensor perturbations would imply a value of Hj high 
enough to make the isocurvature fraction close to unity, which is clearly inconsistent 
with observations [30]. 

Parameters are estimated using standard Markov Chain Monte Carlo techniques, 
as implemented in the publicly available package CosmoMC [31]. As input we use four 
different data sets: 

1. WMAP: the five-year WMAP data [25]. 

2. WMAP+SDSS-LRG+SN: WMAP plus the galaxy power spectrum from the 
luminous red galaxy subsample of the Sloan Digital Sky Survey [26] and the 
luminosity distances of supernovae from the Union compilation [27]. 

3. Planck: Simulated TT, TE and EE spectra up to ^ = 2000 from the Planck 
sateUite [28], assuming isotropic white noise, a sky coverage of 70%, and 14 months 
of observation in the 70, 100, and 143 GHz channels. 

4. CVL: Simulated, noiseless (i.e., cosmic variance hmited) TT, TE and EE spectra 
up to ^ = 2000. 

The mock CMB data sets 3 and 4 are generated using the method discussed in Ref. [29]. 
The hkelihood function C is defined as 

\Ce 



Tr(C7^Q) +ln 



n 



(3.1) 



Xe'ff^-21n£=^(2£+l)/,ky 
e=2 

where Ce = Cf^ with X,Y e [T, E] denotes the mock data covariance matrix, 
Ce = Cf^ + is the total covariance matrix comprising theoretical predictions of 
the CMB anisotropy spectrum Cf"^ and the noise power spectrum N^"^ . The quantity 
n counts the number of observable modes, where n — 2 for observations in temperature 
and £'-type polarisation. 

Using data set 2, we find 95% credible intervals for the observables discussed in 
Sec. 2, 

0.1023 <cuc< 0.1165, (3.2) 
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Figure 3. Degeneracies of ttie isocurvature fraction a witli the other parameters of 
the standard ACDM-model for different data sets. Plotted are the marginalised two- 
dimensional 95% credible contours for WMAP data only (pale red), WMAP+SDSS- 
LRG+SN (red), and for simulated data from Planck (thin pale blue dotted) and a 
cosmic variance limited experiment (thin blue dotted lines). 



3.105 < ln[10^°As] < 3.260, (3.3) 
a < 0.09, (3.4) 

consistent with the findings of Komatsu et al. [20]. As noted by these authors, adding 
large-scale structure and supernova data leads to a significantly tighter bound on a 
because these data sets contain additional information on and thus help alleviate 
the degeneracy problem somewhat. This is illustrated in Fig. 3, where we show the 
two-dimensional joint constraints on a with the other parameters of the vanilla model 
for data sets 1-4. 

For the future experiments and if no isocurvature signal shows up, we forecast 
95%-credible upper limits of 

a < 0.042 (Planck), (3.5) 

a < 0.017 (CVL), (3.6) 

in rough agreement with the CMBpol White Paper [23]. || We emphasise that upcoming 

II The noise properties of the CMB probes considered in Ref. [23] are roughly equivalent to our CVL 
experiment. Some other subtle differences include a slightly larger sky coverage of 80% (which leads to a 
general -^/O.S/O.T— 1 ~ 7% improvement in the parameter sensitivities compared with numbers derived 
from a 70% sky coverage) and the inclusion of _B-polarisation (which has no bearing on isocurvature 
detection). A major difference between Ref. [23] and our analysis is the forecast method: we have 
analysed mock future data, while the forecasts of Ref. [23] are based on the Fisher matrix which 
assumes Gaussian posterior distributions. 
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CMB experiments will break essentially all remaining degeneracies without the need 
for additional data. The improvement one can expect from Planck is partly due to 
this effect, but the improved sensitivity to i?-polarisation will also play a role (Planck 
i?-polarisation data will essentially be cosmic variance limited for I < (9(10)). Going 
from Planck to CVL will increase the sensitivity to a by another factor of 2.5, mostly 
because CVL covers the "sweet spot" in the ^'-polarisation signal-to-noise, where the 
first acoustic peak of the isocurvature signal coincides with the first acoustic trough of 
the adiabatic signal, around I — 200 (see Fig. 2). Because the isocurvature signal affects 
predominantly the largest scales, any significant further increase in the sensitivity to a 
from CMB observations beyond the CVL benchmark value is most likely unreahstic. 

4. Axion pcirameter space 

It is now easy to translate our constraints and sensitivity forecasts on the isocurvature 
fraction a into axion parameters by virtue of equation (2.14). It can be written as 

where we have used our present-day a-constraint from equation (3.2) as a benchmark 
value. Assuming axions are the dark matter, we show this constraint in Fig. 4 with a 
fine marked "Current data." 

In this plot, which is inspired by a similar figure in Ref. [19], we also show the 
relationship between /a and ©i as dashed fines marked with the relevant value of ©i. 
We have also marked the region where the PQ scale is smaller than Hj and where 
therefore our "late inflation scenario" would not apply. Also displayed is the region 
of Hi values excluded by excessive tensor modes, corresponding to r > 0.25. Future 
sensitivities to a from Planck and CVL based on equation (3.5) are shown labelled 
"Planck" and "CVL" respectively 

5. Conclusions 

We have studied current cosmological limits and future sensitivities for axion-type 
isocurvature fluctuations. Our present-day limit on the isocurvature fraction a < 0.09 
(95% C.L.) agrees with similar recent results by other authors. We have used a ACDM 
cosmological model with six standard parameters and in addition a scale-invariant 
isocurvature fluctuation spectrum that is uncorrelated with the dominant adiabatic 
fluctuations. 

The main effect of axion isocurvature fluctuations is to modify the Sachs-Wolfe 
plateau at small CMB multipole orders.^ However, breaking degeneracies with other 
parameters, future CMB probes can signiflcantly improve the sensitivity to a and thus to 

% This phenomenological limitation can be circumvented in a class of extended axion models that lead 
to a significantly blue-tilted spectrum [32], affecting also CMB anisotropics at larger multipoles. 




log(////GeV) 

Figure 4. Exclusion and sensitivity regions in the plane of Hj (Hubble rate during 
inflation) and /a (axion decay constant), assuming axions are all of the dark matter. 
The isocurvature exclusion region based on current data is shown in light blue with 
a boundary given by equation (3.2). The sensitivity forecasts for Planck and CVL 
of equation (3.5) are also indicated. The dashed lines indicate the required Qi for a 
given /a to obtain the full amount of axion dark matter. We also show the region of 
excessive tensor modes and the region /a < Hi where our late-inflation scenario is not 
applicable. 

the presence of axions in the universe. A non-detection of isocurvature modes by Planck 
would improve the limit to a < 0.042, while an ultimate CMB probe limited only by 
cosmic variance in both temperature and E'-polarisation up to £ = 2000 could reach 
a < 0.017, about a factor of five more restrictive than current limits. These forecasts 
agree roughly with previous estimates, notably by the CMBPol Collaboration [23]. The 
small differences are probably due to their using Fisher-matrbc techniques in contrast 
to our full-fledged analysis of mock data sets. 

We emphasise that the CVL sensitivity to a is essentially "as good as it gets", 
since all relevant degeneracies have been lifted and no other cosmological probe besides 
the CMB is directly sensitive to isocurvature fluctuations. In terms of the axion 
decay constant /a, the cosmic-variance limited probe shifts the sensitivity region by 
approximately a factor of five for some fixed Hubble parameter during inflation Hj, and 
thus opens a sliver of parameter space where axions can still leave a detectable imprint 
in the microwave sky. If a tentative signal were found, an experimental search for axion 
dark matter in this parameter range would be crucial, yet technologically challenging. 

The idea of dark-matter axions with /a at the GUT scale or larger has received 
relatively little attention, at least from the perspective of experimental searches, because 
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the initial misalignment angle would have to be much smaller than unity. However, 
in our late-inflation scenario the axion dark matter density is fixed by the random 
number — tt < ©i < +7r with a flat prior distribution. In this case anthropic selection is 
practically unavoidable and implies that the required small Gj value is not unnatural. 
Therefore, as stressed by previous authors, the anthropic axion window is a plausible 
parameter range. In the absence of direct experimental searches in this window, the 
appearance of isocurvature fluctuations is the only opportunity for axions to become 
visible in this scenario. However, in view of the fundmental limitations posed by cosmic 
variance on all cosmological probes, a large chunk of parameter space in /a and Hi will 
likely remain unreachable by this method. 
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